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As is known above 90% of the energy content in Universe is made of unknown dark component. 
Usually this dark fluid is separated into two parts: dark matter and dark energy. However, it may 
be a mixture of these two energy components, or just one exotic unknown fluid. This property 
is dubbed as dark degeneracy. With this motivation, in this paper, a unified dark fluid having 
constant adiabatic sound speed = a, which is in the range [0,1], is studied. At first, via the 
energy conservation equation, its energy density, pd/ Pdo = (1 — Bs) + Bsa~^^^'^'^^ where Bs is 
related to integration constant from energy conservation equation as another model parameter, 
is presented. Then by using Markov Chain Monte Carlo method with currently available cosmic 
observational data sets which include type la supernova Union 2, baryon acoustic oscillation and 
WMAP 7-year data of cosmic background radiation, we show that small values of a are favored in 
this unified dark fluid model. Furthermore, we show that smaller values of q < iO~^ are required 
to match matter (baryon) power spectrum from SDSS DR7. 



I. INTRODUCTION 

Since the discovery of an accelerated expansion of our 
universe [J [3] , a flood of dark energy models have been 
presented to describe an exotic energy component which 
has negative pressure and pushes the universe to an ac- 
celerated expansion. For recent reviews on dark energy, 
please see 0-01 • As is known above 90% of the energy 
content in Universe is made of unknown dark component. 
However, the nature of this dark fluid is still unknown. 
Usually, this dark fluid is divided into two parts: dark 
matter and dark energy. But as pointed out by the au- 
thors [l^, [ll| , there exists another possibility that it is 
a mixture of dark matter and dark energy, or just one 
exotic unknown fluid. This property is dubbed as dark 
degeneracy. A unifled dark fluid is defined from the Ein- 
stein field equation 
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where T°^* is the observed energy components and G^jy 
is the observed geometric structure of Universe. In the 
framework of linear perturbation theory, the micro scale 
properties of a fiuid are characterized by its equation of 
state (EoS) Wd and sound speed c^. As a contrast to the 
reports where a parametrized form of EoS of dark energy 
is assumed in the so-called model-independent formal- 
ism, we study a unified dark fluid with constant adiabatic 
sound speed = a. Actually, the case of zero adiabatic 
sound speed has been studied in Ref. llj, 12;]. And time 
variable sound speed cases were also discussed in [l^ . In 
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this paper, we consider a generalized case^Le cj. = a. 
The similar situation was also discussed in [l4| where an 
effective cosmological constant p\ and dark matter were 
deflned. And it was dubbed as AaCDM model. In fact, 
for a unified dark fluid, one can do different decomposi- 
tion, for example defining dark matter pdm = Pdmoo.~^ 
and the remaining part as dark energy pde ^ Pd ^ Pdm- 
However, without the guide of physics principle, any de- 
composition would be non-proper. Moreover, it may be 
just an entirely whole energy component. And there does 
not exist any decomposition at all. So as a contrast to 
Ref. [l^ , in this paper, we prefer taking it as an entirely 
whole dark fluid and show the properties of dark degen- 
eracy. Furthermore instead of using the location of CMB 
acoustic peaks as done in [l^ . the full information from 
WMAP 7-year data sets will be used in this paper. As a 
result, a tighter constraint will be obtained. 



At first, in section |TT1 the energy density and back- 
ground evolution equation for a unified dark fluid of con- 
stant adiabatic sound speed are shown. The correspond- 
ing scalar flelds action is also presented. In section HIIl by 
using Markov Chain Monte Carlo (MCMC) method with 
currently available cosmic observational data sets which 
include type la supernova Union 2, baryon acoustic oscil- 
lation and WMAP 7- year CMB, we show the parameter 
space. A summary is presented in section ITVl 
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II. GENERALIZED DARK DEGENERACY 
FLUID WITH A CONSTANT ADIABATIC 
SOUND SPEED 



A. Background Equations, Perturbations and 
Instabilities 



Consider a unified dark fluid as a barotropic fluid with 
constant adiabatic sound speed 



= a. 



(2) 



For barotropic fluid the adiabatic sound speed equals the 
speed of which perturbations propagate in the fluid. Gen- 
erally its EoS is given as 



(3) 



where w{p) is the EoS of this unified dark fiuid. From 
the definition of adiabatic sound speed 

dP\ dP dw 
op J ^ dp dp 

after integration, one obtains its solution 



A 

w = a , 

P 



(5) 



where A is the integration constant. The pressure P = 
ap — A can be obtained immediately. Assuming only 
gravitational interaction between energy components and 
considering the conservation of energy of dark fluid 

Pd + iH{l + Wd)pd = 0, (6) 
one has the energy density of dark fluid 
A 



Pd 



l + a 



where B is another integration constant. After the nor- 
malization of scale factor a to 1 today, ao = 1, one has 
dark fluid energy density in the form 



Pd = PdO 



{(1-5.) + 5.0-3(1+")} 



(8) 



where Bs — B/ pdo- In terms of B^ and p^Oj one has 
A = Prfo(l + — Bg). For this unified dark fluid, one 
can do a decomposition, for example, pdm = Pdmoa~^ 
with usual power law and pde — Pd ~ Pdm- If you like, 
you can do many kinds of decomposition. But, without 
any guide of physics principle, any kind of decomposition 
is non-proper. It has the possibility that it is just a whole 
component. And there does not exist any decomposition 
at all. So, in this paper, we will keep it as a whole uni- 
fied dark fluid and do not implement decomposition any 
more. Then one can recast EoS of dark fluid Wd into 

{\ + a){\-B,) 

Wd = a — ^ 57— r. (9) 

(1 - Bs) + S.a-3(i+") ^ ^ ' 

in terms of B^ and a which are model parameters in our 
story. It is interesting that the adiabatic sound speed a 
appears in the power law index of scale factor. It can 
be understood easily as follows. Since = a has re- 
lation with perturbation propagates in fluid, it modifies 
the scaling relation of energy density with scale factor 
naturally. Staring at the expression of energy density of 
this dark fluid, one sees that a cosmological constant is 
recovered when Bs = and dark matter is obtained if 
Bs = 1 and a — are respected. To protect energy den- 
sity from negativity, < i?. < 1 is mandatory. In some 
sense, one can take this unified dark fluid energy density 
as a interpolation from 'dark matter' and cosmological 
constant. Just in this sense, it was dubbed as AaCDM 
model by authors [3]. Then one can find out that the 
EoS of this dark fluid is in the range oi Wd G [— l,a). 
We can show the evolutions of Wd{a) with respect to the 
scale factor a in Figure [U where different values of model 
parameters a and Bg are taken. From the left panel of 
Figure [U one sees that this dark fiuid is dark matter like 
in the early epoch when a approaches to zero in the limit 
(7) a ^ 0. 



J 



Consider a spatially fiat FRW universe 

ds^ = ^dt^ + a^{t) [dr^ + r^{de^ -f sin^ edcj}"^)] , (10) 
I 

H'' = Hi [{I - n^, - ^r) 



where H is the Hubble parameter with its current value 
Hq = lOO/ikm s-^Mpc^i, and fli {i — b,r) are dimen- 



from the Einstein field equations, one has the Friedmann 
equation 

+ Qba-^ + nra-'^Y (11) 

I 

sionless energy parameters of baryon, radiation density 



(l-B3) + B.a-3(i+") 
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FIG. 1. The evolutions of Wd{a) with respect to the scale factor a. In the left panel, Bs = 0.23 is fixed. And the lines from 
the bottom to top correspond to a = 0.01, 0.02, 0.04. In the right panel, a — 0.001 is fixed. The lines from the bottom to top 
correspond to Bs = 0.2, 0.3, 0.4. 



respectively. 

Considering the perturbation in the synchronous gauge 
and the conservation of energy-momentum tensor — 
0, one has the perturbation equations of density contrast 
and velocity divergence for dark fluid 



6d = -(1 + wd)ied + |) - 3H(c2 - wa)6d (12) 



9d = ~n{l - 3cl)ea 



l+Wd 



k^Sd - k^ad (13) 



following the notations of Ma and Bertschinger [l5| . For 
the gauge ready formalism about the perturbation the- 
ory, please see jig . For the dark fluid in this paper, we 
assume the shear perturbation ad = 0. In our calcula- 
tions, the adiabatic initial conditions will be taken. 

The linear and nonlinear instabilities were studied ex- 
tensively in the context of unified dark fluid model [13] ■ 
For the linear perturbations, the adiabatic sound speed 
of the unified dark fluid = a is kept in the range of 
to 1 in our model, so the linear instability is avoided. 
The averaging problem is involved when perturbations 
become nonlinear. The problem coms from the fact that 
(p) ^ pUp)) for unified dark fluid models. For example, 
in the generalized Chaplygin gas (GCG) model [isj 



P' 



(14) 



where A and /3 are positive model parameters in GCG 
model, one can check that 



(p) = -(A//)^-^/(p)^=p((p)), 



(15) 



in the case of /3 7^ 0. However, it is not the issue in 
our model for the linear relations between p and p, i.e. 
p = ap—A. It is a feature of our unified dark fiuid model. 



B. Scalar Fields 

In [l3|, the authors discussed the relations of GCG 
with scalar fields. Then what about the unified dark 
fiuid with a constant adiabatic sound speed? In this sub- 
section, we will derive its Lagrangian density following 
the procedure of 17|. For a sirigle scalar field, the action 
is generalized to the form [l3, [3| 



Sd, — 



(16) 



where £ is the Lagrangian density, (/> is a real scalar field, 
and X is its kinetic term X = — V^V^0/2. Then the 
energy density and pressure are given as 



2Xp^x -P, P = CiX, , 



(17) 



To get the Lagrangian density is to solve the first differ- 
ential equation of Eq. (|17p in terms of X. For our unified 
dark fluid model, one has 



2aXp^x - (1 

For a = 0, the solution is p{X) 
has the solution 

p = aX(i+")/(2")-A/(l + a). 

In terms of X , the energy density is given as 
p = A/(l + a)+X(i+")/(2a)_ 

Thus the Lagrangian density is 

C{X) = p{X) = Q,x(i+")/(2°) - A/(l + a) (21) 

with Q < X < 1/2. One can check that the adiabatic 
sound speed is 



a)p-A = Q. (18) 
A. For a 7^ 0, one 

(19) 
(20) 



P,X 



(22) 



Here please note that it is the adiabatic sound speed not 

vhich i 

filed rest frame, i.e 5<t)\rf = |2Q] . 



the sound speed c^^ = 1 which is defined in the scalar 
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III. CONSTRAINT METHOD AND RESULTS 
A. Implications on CMB anisotropy 

Before going to the issue of cosmic observational con- 
straint, we would like to show some implications on CMB 



anisotropy in this unified dark fluid model. In Figure [2l 
we illustrate how the CMB temperature anisotropics are 
characterized by different values of a, by choosing differ- 
ent values of Bs with the other cosmological parameters 
fixed. 




FIG. 2. The effects on CMB temperature power spectra of a unified dark fluid with constant adiabatic sound speed = a 
characterized by different values of a and Ba- On the left panel, the fixed value of a = is adopted. The bottom red solid 
line is with Bs — 0.45 and top green dotted line is with Bs = 0.35, and the center one is with Bs = 0.40. On the right panel, 
Bs — 0.35 is fixed. The lines take the values of q = O.OOOf , O.OOf , 0.01 from bottom to the top. 



From the left panel of Figure [21 as is expected Bg cor- 
responds to an effective dimensionless energy density of 
dark matter at early epoch. So, increasing the values of 
Eg will make the equality of matter and radiation earlier, 
then the sound horizon is decreased. As a result, the first 
peak is depressed. The variation of the second peak is 
from the variation of ratio of baryon and effective dark 
matter. The right panel of Figure [2] shows the effects 
of a on the CMB power spectra. As shown in Eq. ([5]), 
the values of a describe the possible deviation from stan- 
dard evolution scaling law of effective dark matter. 
It affects the CMB power spectra at large scales I < 100 
via Integrated Sachs- Wolfe (ISW) effect due to the evolu- 
tion of gravitational potential. At small scales, the effect 
comes from the different scaling of effective dark mat- 
ter evolution. From the Figure [H one can read that the 
CMB power spectra are sensitive to the values of a and 
Bs- Then the full CMB information will be helpful to 
constrain the model parameter space. 



B. Method and data points 

To obtain the parameter space, the observational con- 
straints are performed by using Markov Chain Monte 



Carlo (MCMC) method. We modified publicly available 
cosmoMC package |2l|[to include the dark fluid pertur- 
bation in the CAMB [22| code which is used to calculate 
the theoretical CMB power spectrum. The following 7- 
dimensional parameter space is adopted 

P = {LJb, Os, T, a, B„ n,, log[IOi%]} (23) 

where uji, — Q,i,h? is the physical baryon density, Qs 
(multiplied by 100) is the ration of the sound horizon 
and angular diameter distance, r is the optical depth, a 
and Bs are two newly added model parameters related 
to dark fluid, Ug is scalar spectral index. As is the am- 
plitude of of the initial power spectrum. The pivot scale 
of the initial scalar power spectrum kgo — 0.05Mpc^^ is 
used in this paper. The following priors to model pa- 
rameters are adopted: Wf, e [0.005,0.1], 9s G [0.5,10], 
T e [0.01,0.8], a e [0,1], Bs e [0,1], Us e [0.5,1.5], 
log[10"'^°^s] S [2.7,4]. Furthermore, the hard coded prior 
on the comic age fOGyr < < 20Gyr is also imposed. 
Also, the physical baryon density ut, = 0.022 ±0.002 [H 
from big bang nucleosynthesis and new Hubble constant 
Ha = 74.2 ± 3.6kms"^Mpc"^ [H] are adopted. 

To get the distribution of parameters, we calculate the 
total likelihood L oc where is given as 

= XcMB + X^AO + xIn- (24) 
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The 557 Union2 data [281] with systematic errors and 
BAO [25] are used to constrain the background evolu- 
tion, for the detailed description please see Refs. p9l |. 
SN la is used as standard candle. And BAO is used 
as standard ruler. At the last scattering of CMB radia- 
tion, the acoustic oscillation in the baryon-photon fluid 
was frozen and imprinted their signature on the matter 
distribution. The characterized scale of BAO in the ob- 
served galaxy power spectrum is determined by the co- 
moving sound horizon at drag epoch z^i which is shortly 
after photon decoupling. To calculate rs{zii), one needs 
to know the redshift Zd at decoupling epoch and its cor- 
responding sound horizon. In our story, the usual fitting 
formula [26| can not be used for its viability under the 
conditions pb oc a~'^ and pc oc a~^. So, to use the BAO 
information, we obtain the baryon drag epoch redshift Zd 
numerically from the following integration [27j 



dri'id 



da R 



1 



(25) 



where R = 3pb/4:pj, ut is the Thomson cross-section and 
Xe{z) is the fraction of free electrons. Then the sound 
horizon is 



0.01, the global fitting results are summarized in Ta- 
ble H] and Figure [31 We find that the minimum is 
Xmm = 8009.424 which is a little bigger than that for 
ACDM model 8009.116 for the same cosmic observational 
data sets combination. From the Table HI one can see 
that small values of = a are favored in this unified 
dark fluid model when the WMAP 7-year observation, 
SN Union2 and BAO data points are used as cosmic con- 
straint. It is worth to mention that our result is consis- 
tent with that reported in Ref. [3] where a relative loose 
constraint a = 0.01 ± 0.02 was given. Correspondingly, 
we plot the CMB power spectrum for the mean values 
estimated from MCMC analysis in Figure IH where the 
observational data points of 7- year WMAP with uncer- 
tainties are also included. As comparisons, the CMB 
temperature power spectrum for ACDM with same cos- 
mic observational data sets combination and results from 
[30j are also shown. One can see that the CMB power 
spectrum for this unified dark fluid model is well inside 
the error bars of the binned measurements from WMAP 
7- year results and almost matches to ACDM model. This 
strongly implies that current cosmic observational data 
combinations of SN Union2, BAO and WMAP 7- year can 
not almost discriminate a unified dark fluid with con- 
stant adiabatic sound speed from ACDM model, though 
ACDM model is still slightly favored. 



rs{zd) 



drjCsil + z). 



(26) 







where Cs = l/\/3(l -I- R) is the sound speed. Also, to 
obtain unbiased parameter and error estimates, we use 
the substitution 123 



dr—, — -rs\Zd), 



fs{zd) 



(27) 



where dz = rs(zd)lDy{z)., fs is evaluated for the fiducial 
cosmology of Ref. [2^ , and Zd is redshift of drag epoch 
obtained by using the fitting formula [2^ for the fidu- 
cial cosmology. Here Dy(z) = [(1 + z)^D^c2:/i?(z)]^/^ is 
the 'volume distance' with the angular diameter distance 
Da- In this paper, for BAO information, the SDSS data 
points from [25 are used. For CMB data set, the tem- 
perature power spectrum from WMAP 7-year data [SO^ 
are employed as dynamic constraint. 



Prameters 



Mean with errors 



0.0225 



+0.000504+0.000979 



- +0.00246+0.00491 



0.0874 



+0.00607+0.0243 
0.007010.0225 



log[10i°A] 

Age/Gyr 

Ho 



n nnn/ia7+0-OOOii'i'+o. 000728 

U.UUUy:0 / -0.000487-0.000487 
n OOQ+0. 0134+0. 0274 
"•^^»_0. 0133-0. 0249 
n nyc+0. 0132+0. 0265 
U.y/O_o. 0131-0. 0249 
n |-,oi +0.0331+0.0663 
O.U01_Q 0327-0.0616 
n QKfi+0. 00152+0. 00298 
U.yOD_o 00150-0.00300 
1 fiQQ+0. 108+0. 210 
iO.Oyy_o. 1085-0. 215 



044+°-''''-'^''+''-''''^'''' 
10.523"* 



71.341 



hi. 1948+2. 393 
-1.1618-2.299 
+ 1.278+2.573 
-1.268-2.478 



TABLE I. The mean values of model parameters with la and 
2a errors, where WMAP 7-year, SN Union2 and BAO data 
sets are used. 



C. Fitting Results and discussion 



After running 8 independent chains and checking the 
convergence to stop sampling when i? — 1 is of order 
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FIG. 3. The ID marginalized distribution on individual parameters and 2D contours with 68% C.L. and 95% C.L. by using 
CMB+BAO+SN data points. The shade regions show the mean likelihood of the samples. 



Up to now, by using the geometric informations from 
SN, BAO and CMB data sets, one has obtained a rela- 
tive tight constraint to the model parameter space. How- 
ever, for any real Universe model, one needs to consider 
the large scale structure formation. That is to say, the 
dynamic evolution information from galaxy formation 
would be employed. It is not an easy issue for unified 



dark fluid model. Here, we just show that the smaller 
values of a < 10~^ are required for matching the mat- 
ter (baryon) power spectra data points from SDSS DR7 
[31} , please see Figure [5] It means that when the dy- 
namic information is used, for example large scale struc- 
ture formation, a tighter constraint will be obtained. For 
the constraint to this unified dark fluid model from large 
structure information, we leave it for the future work. 



IV. SUMMARY 



In this paper, we have studied a unified dark fluid 
model with constant adiabatic sound speed = a. This 
property is called dark degeneracy. By using MCMC 
method with current available cosmic observational data 
sets combinations of SN Unioin2, BAO and WMAP 7- 



year results, we find that small values of a are favored 
in this unified dark fluid model as shown in Table U and 
Figure [31 For small values of a, this unified dark fluid 
model approaches to ACDM model closely. Our analysis 
shows that ACDM model is still slightly favored when 
SN Union2, BAO and WMAP 7-year observational re- 
sults as useful constraint. But, there exists the possibil- 
ity that the unknown above 90% content of our universe 



7 




FIG. 4. The CMB C; power spectrum v.s. multiple moment I, where the black dots with error bars denote the observed 
data with their corresponding uncertainties from WMAP 7-year results, the red solid line is for the unified dark fluid model 
with mean values as shown in Table [H the blue dashed line is for ACDM model with mean values for the same data points 
combination. And the green doted line is for ACDM model with mean values taken from with WMAP+BAO+Ho constraint 
results. 



is made of a unified dark fluid, the so-called dark degen- 
eracy property. To break the degeneracy, more cosmic 
observational data sets would be included. We also show 
that smaller values of a < 10~^ are required to match 
the matter (baryon) power spectra from SDSS DR7. We 
expect our study can shed light on the understanding the 
dark side of our Universe. 
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